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ABSTRACT 

We present the density structure of dark matter in the outer parts (to about 7 effective radii) of 
the eUiptical galaxy NGC 4636, from the radial velocities of 174 globular clusters, using the non- 
parametric, inverse algorithm CHASSIS. We find the galaxy to be rich in dark matter, with i?-band 
mass-to-light (M/L) ratios rising to about 30, at nearly 4i?e; the if-band M/L at about 3i?e is found 
to be nearly 10. The result does not depend on applying the method to the red and blue globular 
clusters separately. This estimate of M/L is higher than the previous analysis from the same kinematic 
data. We also find that the dark matter distribution is highly concentrated towards the inner halo. 

Subject headings: galaxy: kinematics and dynamics — galaxy: globular clusters: individual (NGC 4636) 



1. INTRODUCTION 

The overwhelming presence of dark matter in the out- 
skirts of disk galaxies has been invoked to explain flat ro- 
tation curves out to se veral t imes the optical d imensions 
of the se systems, (e.g. . iRoberts. 1969, : iFaber fc Gallagher 
Il979| ). However, the inferred dark matter content of 
cores of galaxies, in many cases, seems to fall well short 
of the standard model p redictions (e.g., iGentile et al.l 
120041 : iFerreras et al.l[2007f ). In the same spirit, intrigu- 
ing studies of dark haloes of early-type galaxies, using 
planetary nebulae as tracers of the distribution of mat- 
ter, seem to indicate the under-ab undance of dark mat- 
ter (e.g.. IRomanowskv et al.l[2003f l. This view has been 
challenged, however, on grounds that the subtleties in 
the very processi ng of the kinematic data can lead to 
spurious answe rs jPekel et al. 120051: 1 Sambhus et al.ll200l 
iDouglas et al.l [20071 ). In particular, the mass-anisotropy 
degeneracy is a difficult issue to resolve, given the limited 
breadth of available kinematic measurements. 

In spite of this shortcoming, it is highly appealing 
to use kinematic information in verifying that the dark 
matter distribution of elliptical galaxies is in concor- 
dance with with the predictions of cosmological simu- 
lations, performed within the ACDM paradigm. Glob- 
ular clusters provide an independent source of test 
particles, since their formation histories are differ- 
ent from those of planetary nebulae, and the nature 
and distribution of their orbits would provide inde- 
pendent information about the distribution of the un- 
derlying dark matter. Mass modeling based on glob- 
ular cluster dyna mics has been tried out on nearb:v 



early-type galaxies (*C6 te et al]|2003t 



Schubcrth et al. 2006; Woodlev et al 



Bridges et al 



r ear by 

Mm 



120071 : iHwang et all 



20081: IRichtler et alM206§ ). and in many cases there have 



been appreciable discrepancies between these results and 
the models obtained from planetary nebulae. Such en- 
deavors highlight the necessity of not only studying a 
large number of systems for both sets of tracers, but also 
investigating the dependence on models, as well as trying 
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out different methods of mass reconstruction, in partic- 
ular those that do not explicitly depend on the charac- 
terization of models. 

NGC 4636 is an example of an early-type galaxy for 
which the different modes of mass modeling have led to 
a large discrepancy in the recovered mass, even at the ef- 
fective radius (i?e=101".7, Schuberth et. al 2006). NGC 
4636 lies near the southern border of the Virgo cluster 
and is relatively isolated. It has a radial velocity similar 
to that of the Virgo cluster, even though it is about 3 Mpc 
from the center of the cluster. Furtherm ore, the galaxy 
lies at the center of a poor g roup (Osmond fc PonmanI 
120041: iMiles et al.ll2004 |2006[ ). possibly falHng into the 
cluster. Its unusual properties have attracted detailed 
multi-wavelength research for several decades. For in- 
stance, NGC 4636 is found to be very bright in X-rays 
(Lx ~ 10'*^ergs/s), with unusual features in the ho t ISM 



dForman et al. 1985; Matsushita ct al. 1998; Jones et al 
'2002; Loewenstein fc Mushotzkv. ,2003 : ,0'Sullivan et al 
2005). 

NGC 4636 has an anomalously large specific frequency 
of globular clusters (jPirs ch ct al. 2005;), comparable to 
that of central galaxie s of far richer syste m s, like that of 
NGC 1399 in Fornax (jDirsch et al.ll2OO30 . IDirsch et all 
(|2005l ) provide details of several thousand globular clus- 
ter candidates around this galax y, of which rad ial veloci- 
ties of 174 GCs were presented in lSchuberth et a l. (2006), 
where a dynamical analysis of this data, on the basis of 
Jeans Equation, was also performed, for assumed values 
of the stellar mass to light ratio {M/L) and anisotropy. 
Naturally, the recovered mass profile is the projected one, 
which only provides a lower limit to the actual mass dis- 
tribution. In this approach, the total mass density profile 
needs to be parameterized (as an NEW model, say) and 
constraints are recovered for these characteristic parame- 
ters. Thus, this method is not well-suited for the recovery 
of the distribution of total mass. Further, the estimation 
of the density parameters is awkwardly sensitive to the 
choice of the stellar M /L that is used to scale the lumi- 
nosity density prof ile: ISchuberth et al.l (I2006D use M/L 
values obt ained by [B ell et al.' ('2003'), iKrona witter et al.l 
(pOOO) and iLoewens tcin & Mushotzkv (2003), from var- 
ious samples of elliptical galaxies. There is of course no 
reason why this ratio has to be a constant and not vary 
with radius in the radial range covered by the kinematic 
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data (projected r adius R <30 kpc ). 

According to ISchuberth at al.1 ()2006l ). the range 
of mass profiles of the halo of NGC 4636, corre- 
spond ing to an assumed NFW density (jNavarro et al.l 
19961) ■ do indeed straddle the distribution indicated by 



Matsushita et alJ (fl998,'). though the profile obtained by 



Loewenstein fc Mushotzkvl (I2003D is f ound t o be too mas- 
sive. even for the ISchuberth et all (|2006f ) model that 
corresponds to the highest dark matter content within 
30 kpc. 

Clearly, the use of different methods, with varying 
model assumptions, on different data sets, disagree in 
their attempts to find the underlying dark matter distri- 
bution within the same volume. Here, we take an ap- 
proach which does not assume a parameterized model. 
We u se the algorithm CHASSIS (jChakrabartv fc Sahal 
1200 It) to analyze the available kinematics of the same 
sample of 17 4 GCs that i s used in the Jeans Equation 
approach of iDirsch et al.l (|2005f) . This algorithm has 
been calibrated agains t the N-body real ization of two 
star clusters ([Chakrabartv fc Portegies Zwart 200 41 , and 
applied to estimate the central ma ss structure of the 
Galaxv (iChakrabartv fc Sahal [2f)0lh and the GC M15 
(|Chakrabartvll2006H . It is also used to estimate the dark 
matter content in NGC 3379, using planetary nebulae 
kinematics (Chakrabarty 2008, in preparation). 

The basic formalism of CHASSIS is discussed in the 
following section. Section 3 describes the observational 
details of the GCs used in this analysis, and in Section 4, 
we present the results obtained from CHASSIS for the 
AI/L distribution in the outer parts of NGC 4636, and 
compare our results with previous work. Section 5 sum- 
marizes the implications of our results in the determina- 
tion of the content of dark matter in the outer halo of 
NGC 4636. 

2. THE CHASSIS ALGORITHM 

CHASSIS works under the assumption that the input 
kinematic data is drawn from an equilibrium phase space 
distribution function that is isotropic in nature. Also, 
the system geometry is assumed to be spherical. The 
algorithm produces a pair of functions: the equilibrium 
phase space density and the (total) mass density, that 
best describe the observed data. These two character- 
istic functions are sought simultaneously, using a max- 
imum likelihood approach that employs a sophisticated 
optimizer - the Metropolis algorithm. The inputs to the 
algorithm are the position on the plane of the sky, and at 
least one velocity component (usually the radial velocity 
Vz) of individual GCs in the system. 

As in all recursive algorithms, CHASSIS too requires 
initial guesses (or seeds) for the answers it seeks, namely 
the distribution function and mass density profile (repre- 
sented in the code as histograms over energy and radius, 
respectively). The final answer should be insensitive to 
the choice of this guess; robustness checks are carried out 
to confirm this with the kinematic data used. 

During any run, at the end of the first step, the seeds 
for the mass density and the distribution function are 
slightly tweaked. Likewise, at the end of subsequent step, 
the profiles are modified, (both in shape and overall am- 
plitude), over their previous forms, subject to the con- 
straints of monotonicity and positivity. This is carried 
on until the global maximum in the likelihood function 



is identified. 

The ±1-CT spread in the sample of the density and dis- 
tribution functions, that correspond to the neighborhood 
of this global maximum, readily provides the 1-a error 
bars on the recovered profiles. Quantities that are esti- 
mated from the recovered profiles, such as the enclosed 
mass profile and the velocity dispersion profile, bear the 
signature of this extent of error. It may be noted that 
these errors stem from the uncertainties in identifying 
the global maximum in the likelihood function, and are 
essentially different from the observational errors. The 
errors in the velocity measurements are incorporated 
into the analysis by convolving the projected distribu- 
tion function with the distribution of the observational 
errors, (assumed Gaussia n). Further det ails about this 
algorithm can be found in IChakrabartv fc Saha (2001]). 

CHASSIS works by projecting the distribution func- 
tion in each step, at the current choice of the potential 
(calculated from the current choice of the mass density) , 
into the space of observables. The product of the result- 
ing projected distribution functions, corresponding to all 
the (rp, Vz) pairs in the data set, defines the likelihood 
function. 

3. DATA 

The positions and radial velocities (vz) for the 174 
globular c lusters that ar e used in this study, (adopted 
from iDirsch et all (|2005f) ). are sho wn in Fig. [J Based 
on surface brightness fluctuations, iTonrv et al.l ()200lD 
quotes a d istance to NGC 4636 to be 14.7 Mpc, while 
IDirsch e~al . (2005) quote a distance of 17.7 Mpc based 
on the peak of the GC luminosity function. Both mea- 
sures have considerable uncertainty, so we adopt a dis- 
tance of 16 Mpc to NCG 4636, which translates to 
1" =77.5 pc, approximately. 

As has now been found in a wide range of early-type 
galaxies, the colors of the GCs form a bimodal dis- 
tribution, with the redder, more metal-rich G Cs more 
likely to harbor low-mass X-ray binaries (e.g. [jprdani 
2004; 'Posson-BrowneLaU l200a iKundu fc ZepJ 120071 : 
Woodlcy et al. 200^, there being significant differences 
in key structural as well as chernical p roperties in the two 
populations f e.g. Jordan et al.ll2007f) . According to the 
photometry of lDirsch et al.l ( 2005( ). the blue and red GCs 
of NGC 4636 can be characterized as C — i? color being 
greater than or less than 1.55 respectively. We show this 
classification in the left panel of Figure |T] Another way 
of classifying the data at hand is along the lines of the 
magnitude of the measured velocity errors; any GC with 
a measurement error exceeding 35 kms~^ is assigned to 
one kinematic class while those with higher velocity er- 
rors for the other group. We carry out runs with all 
GCs, and separately with the sub-samples that charac- 
terize each of these photometric and kinematic classes. 

The right panel of Fig. [T] shows the spatial distribu- 
tion of the globular clusters used in this study. This 



is a small fraction of the GC candidates of Dirsch et al 



20051) . and represents the clusters that iSchuberth et al 



2006D report as belonging to the galaxy, on the basis 



of spectroscopic data. As is apparent from this picture, 
the number of observed GCs decreases rapidly beyond 
about 30 kpc, and beyond 50 kpc, there are only 2 ve- 
locity measurements available. This radial profile of the 
measured data set is important in determining the choice 
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Fig. 1. — Left: Measured radial velocities of the globular clusters of NGC 4636 used in this analysis, as a function of projected distance 
from the center of the galaxy. The red points represent the 76 red GCs {C — R color > 1.55), while the blue points the 98 blue GCs 
(C — i?< 1.55). Right: The distribution of these GCs (blue crosses) on the sky, superposed on a optical DSS image of NGC 4636. The two 
yellow circles correspond to radii of 15 and 30 kpc respectively. 



of the radial binning adopted in this analysis. 

For the application of CHASSIS, the radial bin width 
needs to be such that the left edge of the innermost bin 
must not exceed the smallest projected radius at which 
the radial velocity data is available. Working with too 
small bin widths would lead to bins that contain no veloc- 
ity information at small radii. On the contrary, adopting 
too large a bin width can lead to spurious density pro- 
files, especially nearer to the center of the system, where 
the gradient in density is higher than in the outer parts. 
Given these constraints, we experimented with bin sizes 
and found that a choice of 2 kpc is adequate for this 
data set. This value allows us to span radial distances 
over a range of 4.2 kpc to about 55 kpc. However, the 
results recovered by CHASSIS are not crucially sensitive 
to the exact choice of the bin width; this feature of the 
code will be demonstrated below with results from ex- 
periments done with assorted bin width values. 

The 3-D total mass density that is recovered by CHAS- 
SIS is then used to calculate the enclosed mass profile 
which is subsequently compared to the cumulative light 
profile, in order to extract the mass-to-light ratio {M/L) 
profile. In this context, we use the fT-band surface 
brightness profile of NGC 4636, obtained f rom the pho- 
tometry of the 2MASS large galaxy Atlas (jJarrett et al.l 
[2003) . kindly provided to us by Tom Jarrett. This is 
shown in red in the left panel of Fig. [2]). We deproject 
this under the assumption of spherical symmetry (right 
panel of Fi g.[2l) with a non-parametric dep rojection code 
DOPING (ChakrabartxiLFerraesc '2008^ We a lso use 
the i?-band photometry from Dirsch et all (|2005f) to es- 
timate the M/L profile in the i?-band. 

4. RESULTS 

To begin with, the inscnsitivity of the results to the 
choice of the first guess for the distribution function 
needs to be established. For this purpose, we undertake 
to parameterize the seed for the density by an NFW-like 



profile (jNavarro et al.lll996D : 

(r/rc)"(l + r/rcY 

The phase space distribution function is either held as 
a power-law (with a power-law index of /3) or as an ex- 
ponential of the effective energy (e). Thus, the seed for 
the density distribution is characterized by a total of 3 
parameters, namely po, rc, a, out of which, the ampli- 
tude or the central density parameter po is found to have 
no effect whatsoever on the result. Three different seeds 
characterized by the form of the DF and values of Tc, a, 
were implemented in three distinct runs: 

• RUN I - the seed distribution function is set to 
exp(e), rc—50 kpc and a=2.8; kinematic data of 
all 174 GCs are used. 

• RUN II - the seed distribution function cx e^, rc=5 
kpc and Q!=1.8; all GCs are used. 

• RUN HI " the seed distribution function is oc e, 
rc=15 kpc and a=2.3; again, the whole sample is 
used. 

The density profiles resulting from these three runs are 
shown in the lower left panel of Fig. [3l As indicated by 
the results from the different runs, the recovered density 
profiles overlap within the error bars, at almost all radii. 
This lends confidence in the functionality of CHASSIS. 
The estimated density profiles are then used to calculate 
the enclosed mass distributions from which (lower right 
in Fig. [3]), the radial distributions of M/L are recovered 
in the R and K bands (upper left Fig. [3]). Since the 
i^-band photometry is available to about 23.5 kpc, the 
run of M/L is also displayed till this radius, while the 
i?-band M/L is limited to 30 kpc. As expected from 
a typical color of an early-type galaxy, the lower limit 
on the ratio between the i?-band and X-band M/ L val- 
ues, at 20 kpc, is about 2.5. The DF recovered from 
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Fig. 2. — The /?-band (lower curve; from Dirsch et. al 2005) and _R'-band (upper curve; 2MASS. [Jarrett et al.l (200^ ) surface brightness 
profiles of NGC 4636 along the major axis, are shown in red on the left. The deprojected luminosity density profiles (upper curve for 
ii'-band and lower curve for i?-band) are shown on the right and the projections of the same have been overlaid in black over the brightness 
profiles, in the left panel. The deprojection was carried out assuming that the intrinsic geometry of the tracer population is spherical. 



three of the runs is used to calculate the velocity disper- 
sion along the line-of-sight. These have been represented 
in the upper right in Fig. [3] (It is worth pointing out 
here that in t he analysis of Chandra and XMM-Newton 
observations (|0'Sullivan et al.ll2005f ). it has been found 
that between 25-30 kpc, there is an abrupt transition 
between the hot interstellar medium of the galaxy to the 
intergalactic medium of the surrounding group). The ro- 
bustness of the algorithm to the choice of bin width is 
displayed in Figure ID 

To confirm that the results of RUN I, RUN II and 
RUN III are not artifacts of the details of the sampling of 
the tracer population, we performed individual runs with 
sub-samples that correspond to different GC color and 
(radial velocity) measurement accuracy classes. These 
runs are described below: 

• RUN IV - the seed distribution function is set to 
the same as for RUN I; kinematic data of the 76 
red GCs are used. 

• RUN V - the seed is as in RUN I; the 98 blue GCs 
are used. 

• RUN VI - again, the seed is identical to that used in 
RUN I; implemented velocity data are that of the 
121 GCs (red and blue) that hav e velocity errors 
lower than a cutoff of 35 kms~^ (jSchuberth et al.l 
[2mfil . 

These three runs are carried out with bin widths of 2 kpc. 
The size of the error bars on the recovered DF and den- 
sity distributions will be higher for these runs than for 
RUNS I, II and III, which employ the whole sample, i.e. 
a greater number of data points. The results are shown 
in Fig. m 

It is notable that the division of the whole sample of 
the GCs by color did not yield significantly different total 
mass densities; there is indeed a trend for the bluer GCs 
in the sample to be on the higher side of the mass of the 
redder GCs, but as is apparent from Fig. [5l the distinc- 
tion in the recovered mass profile is not strong enough 



to be deemed significant, at ±l-a level. The effect of 
division by kinematic accuracy is even less potent. 

Thus, it appears that the M / Lk distribution within 
the first three effective radii, can at most be about 12 in 
the iiT-band, when the errors of the deprojection, as per- 
formed by DOPING are ignored; when the errors of the 
deprojection are included, the range of the M/L value 
recovered at about 23.5 kpc (« Rf.) is found to widen to 
about 8 to 14. 

4.1. Effect of Assumptions on Mass Recovery 

Within CHASSIS, the effect of imposing isotropy in ve- 
locity space, where anisotropy prev ails, is to recover spu - 
riously high mass density values (jChakrabartvl l2006f l. 
Though in general it is difficult to quantify the extent 
of this bias since it depends on the structure of the un- 
known distribution function, in the case of NGC 4636, 
we can surmise that the effect of anisotropy is not signif- 
icant at th e 1-a level. We i n fer th is using the following 
result that ISchuberth et al.l (|2006( ) report: it is the red 
sub-sample of the observed GCs that were found to ex- 
hibit "significant rotation" in contrast to the blue GCs 
which do not show a "significant signal" . However, from 
Figure [5] we see that the mass density profiles that we 
recover for the red and blue sub-samples are consistent 
with each other within the error bars (though the density 
distribution of the red GCs is on the higher side). 

As for the d e viatio n from sphericity is concerned, 
ISchuberth et al.l (|2006f) suggest an ellipticity of 0.15 for 
NGC 4636. This is in reference to the photometric ap- 
pearance of the galaxy while observational constraints 
on the spatial distribution of the dark matter halo of 
the galaxy, within which the gl obular cluster sys t em re - 
sides, is harder to come by3. ISchuberth et~an (l2006h 
treat NGC 4636 as spherical, motivated by the confirma- 
tion of "modest" deviations from sphericity in the inner 
parts of the galaxy (courtesy Dirsch et. al 2005), where 
most of the GCs live. 

However, it needs to be appreciated that in CHASSIS, 
the assumption of sphericity pertains to the geometry 
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Fig. 3. — Lower left: The recovered 3-D total (dark+luminous) density profiles of NGC 4636, to about 7 times the effective radius (which 
is 101". 7 ~ 7.8 kpc), from three different runs, that were performed with distinct forms of the initial guesses (seeds) for the DF and mass 
density distributions. Results for RUN I are in black, RUN II in red and RUN III in green (details in §4). Lower right: Enclosed mass 
profiles from the density distributions recovered from these runs. The dotted line in black represents 30 kpc, the radius at which the mass 
predicted by CHASSIS is compared to the mass estimates from other work. The dashed line represents the 5i?,e mark. Upper left: The 
ii'-band (lower curve) and i?-band (upper curve) mass-to-light ratios in the outer parts of NGC 4636. Note that the M/L distributions 
in this figure do not include the uncertainties involved in the deprojection of the observed brightness profile into the intrinsic luminosity 
profile, in order to represent the extent of the errors in the analysis due to CHASSIS alone; allowing for these errors would widen the 
range of T k at about 24 kpc (edge of the available photometry) from about 9-12 to about 8-14. The quality of the ij-band photometry 
implied truncating the M/Ln estimate to within 30 kpc even though mass estimates are available to larger distances. Upper right: The 
velocity dispersion along the line-of-sight, estimated from the isotropic equilibrium distribution function recovered by CHASSIS, from the 
runs RUN I, RUN II and RUN III. 



of the total gravitational potential that the set of test 
particles sit in, rather than to the spatial distribution 
of the test particles. Thus, we assume that it is the 
outer parts of NGC 4636 (including its dark halo) that is 
spherical. In any case, mistaking an ellipsoidal system as 
spherical will imply overestimation of the enclosed mass 
profile that is calculated from the recovered mas density 
distribution. Erroneous formulation of the potential will 
also add uncertainty to the recovery of the mass density 
distribution itself though seeking to quantify this bias 
is again not direct, given that it depends on the very 
unknown that we are trying to constrain, namely the 



mass distribution in this system. 

4.2. Comparison with Other Work 

The hot X-ray emitting in terstellar medium 
(ISM) has been well-s tu died (iJones et all 120021: 
'Loew enstein fc Mushot"zkvl 120031 : lO'Sullivan et all 
2005;), with the deepest observations going out to about 
30 kpc which is a radial extent that is well sampled 
by the GCs used in this analysis. [Matsushita et al.l 

^ The globular cluster system in M87 was reporte d to be el- 
liptica l in projection (McLaughlin et al. 1994), though [Cote et al.l 
l|2001h approximate the 3-D spatial distribution of this system as 
spherical. 
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Fig. 4. — This shows the concordance (within ztl-cr error bars) of 
density profiles recovered from runs performed with initial condi- 
tions corresponding to RUN I and bin widths of 2 kpc (in black), 1 
kpc (in red) and 5 kpc (in green). The unavailability of kinematic 
information in 1 kpc sized bins, beyond about 30 kpc, limits the 
radial range over which the density can be sought in this case. 

(|1998D and lLoewenstein fc Mushotzkvl (|2003[ ) concluded 
that NGC 4636 is an extremely dark matter dominated 
galaxy, on the basis of their analysis of ASCA and 
Chandra X-ray observations respectively, by fitting hot 
plasma models and assuming hydrostatic equilibrium. 
When compared carefully, their mass profiles appear 
to be significantly di fferent from each othe r - the 
enclosed mass profile of [Matsushita et all ()1998l ) flattens 
out beyo nd 10 kpc, reaching 5 x 10|^ Mp, at 30 Kpc, 
while the iLoewenstein fc Mushotzkvl (|2003f ) profile rises 
steadily to beyond 10^^ Mq at the same distance. 

However, as is seen from the detailed an alysis of the 
Chand ra a nd XMM-Newton obser vations bv lJones et al.l 
(|2002D and lO'Sullivan et al.l (|2005f ). the hot gas has large- 
scale arm like features that are indicative of AGN out- 
burst, and thus the equilibrium assumption might not 
be a good one. These observations also indicate that 
at about 30 kpc, there is an abrupt transition of the hot 
interstellar medium of the galaxy NGC 4636 to the intra- 
cluster medium of the group surrounding it, where the 
gas is of significantly lower temperature and abundance. 
From these observations, it is hard to decouple the dark 
matter halo of the surrounding group from the dark halo 
of the galaxy itself, but inner to 30 kpc, the dark halo of 
the galaxy must dominate. 

In terms of the shape of the predic ted mass profiles, 
we fin d that unlike the mass profiles of iMatsushita et "all 
(|1998l ). (which exhibit a distinct flattening outside an 
intermediate radius of ~10 kpc), our cumulative mass 
profiles do not indicate any sharp flattening - instead, 
the slope of the profile gradually tapers off from about 
10 kpc, such that even at about 50 kpc, it is only slightly 
rising. 

This is similar in natur e to the mass profile obtained 
bv lSchuberth et al.l ()2006f ) from their Model 14, i.e. their 
mod el with the highe st enclosed mass within 30 kpc. In 
fact. lSchuberth et al.l(l2006il specula te that the flattening 
observed bv lMatsushita et al.l (|1998f ) could be a pecuhar- 



ity of the ASCA data itself. In any case, the apparent 
lack of hydrostatic equilibrium indicates that X-ray de- 
rived mass profiles might not be reliable. This may ex- 
plain why our estimate of the mass of NGC 4 636 falls 
short of what Locw enstein fc Mushotzkvl ()2003f ) suggest 
on the basis of their X-ray studies. 

The mass profi les inferred from a Jea ns analysis of the 
GC kinematics (jSchuberth et al.l |2006[ ) also yield high 
values of the dark matter content in this region - the 
highest and lowest values from their models yield a cumu- 
lative mass of 0.5-0.95x10^^ M© at 30 kpc. At the error 
level of +l-(7, our recovered ma ss at about 30 kpc is in 
excess of that suggested by the ISchuberth et al.l (|2006l ) 
analysis, from the same observational data, the lowest 
range of our errors being just about consistent with the 
highest range obtained by them. Here we recall that 
over-estimation of mass density, at radii where velocity 
anisotropy is present, is expected to be an artifact of 
the assumption of isotropy (though as discussed in Sec- 
tion 4.1, the effect of this artifact is not expected to be 
significant). Thus, the mass-to- light ratios presented in 
this paper are expected to be on the higher side of the 
true values in the two bands discussed {K and R). The 
recovered velocity disper sion profile embr a ces th e disper- 
sion values depicted bv ISchuberth et all (120061 ). within 
the ±l-a error bars. 

It needs to be appreciated that our work is not reliant 
on photometry for the extraction of the all important 
total mass density of the system, and thus we are able to 
offer mass estimates much further into the halo of NGC 
4636, than any of the earlier attempts. We can, in fact, 
calculate the mass density to just inside 7 Re- Of course, 
the computation of the M/L profiles is constrained by 
the limitations of photometry. 

5. DISCUSSION: DARK MATTER IN THE HALO 
OF NGC 4636 

We compute the circular velocity profile of the GC sys- 
tem of NGC 4636 to understand the distribution of dark 
matter in this system. This Vc profile is shown in black, 
as a function of radius in the left panel of Figure [6l The 
profile in red is given as the function a/ GM^^^/r, where 
G is the universal gravitational constant and Msr^ is the 
mass that is enclosed within 3 effective radii. Thus, the 
function plotted in red merely exhibits a Keplerian r^^/^ 
fall-off with radius. The errors on this function are the 
±1 — cr errors on Msr^, as recovered by CHASSIS. 

As is apparent from this figure, there is no significant 
difference between the black and red profiles to about 45 
kpc. This indicates that the dark halo contribution to 
the mass, in the radial range 3i?e to about 5.8i?e is not 
significantly different from the mass at 3i?e. This lower 
radius would be normally expected to bear a lower dark 
matter contribution than the part of the galaxy outside 
5i?e- That this expected trend is not statistically valid 
for the case of NGC 4636, can be interpreted to conclude 
that the distribution of mass in the dark halo of this 
galaxy is highly concentrated. This explanation suggests 
itself readily, since the prospect of low dark mass content 
can be ruled out on the basis of the very high AI / L ratio 
values that we have recovered. 

Driven by this hu nch, we proceeded to fit an NFW 
density distribution (jNavarro et al.l 0*9961 ) to the density 
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Fig. 5. — Top Panel: Total mass density, enclosed mass and M/L in the R (upper curve) and if-band (lower curve), from kinematic 
data of all clusters (RUN I, in black), red GCs (RUN IV, in pink) and blue GCs (RUN V, in blue). Bottom Panel: Total mass density, 
enclosed mass and M/L in the R (upper curve) and /i'-bands (lower curve) from kinematic data of all clusters (RUN I, in black) and the 
GCs with the relatively better kinematic accuracy (RUN VI, in red). 



that is recovered by CHASSIS, from one of the presented 
runs (say, RUN I). Of course, this fit is non-trivial and 
depends on the radial range in the data over which the 
fit is sought. We performed a recursive routine to con- 
strain the mass {Mg) and length scales {ts) in the NFW 
form, such that the extracted rg lay within the radial 
range under consideration. The iterations were carried 
out till convergence was spotted within the errors of the 
fit, which were high. The radial range over which the fit 
was sought extended from 4i?e to about 7i?e while the 
values of and Mg were recovered to be 33.7 kpc±ll% 
and about 1.7x10^^ Mq ± 8%, respectively. This fit is 
shown, superimposed in red, on the recovered density (in 
black) , in the right panel of Figure O This recovered 
would indicate a very high concentration of at least 9, 
for an assumed M200 of at least 2x10^^ M©. This is 
i ndicative of a s t rongl y concentrated system. 

iLintott et al.l (|2006f ) have found that in a sample of 
about 2000 ellipticals from the Sloan Digital Sky Survey, 
the distribution of c is log-normal, with the peak lying 
in the range of 3-10. Thus, on this scale, the dark halo 
of NGC 4636 is predicted to be very highly concentrated 
indeed! 

Recently, iNapolitano et al.l (|2007( ) have suggested a di- 



chotomy in the dark matter distribution of nearby early- 
type galaxies (probed by the Planetary Nebula Spectro- 
graph) , where one class (the "ordinary" , fast-rotating, 
discy/cuspy, early type systems) exhibit rapidly drop- 
ping velocity curves while the other slowly rotating, 
boxy/cored systems display flatter trends in their rota- 
tion curves. From Figure [5] we can see that the medial 
value of the circular velocity increases from about 400 
km s~^ to an intermediate peak of about 450 km s^^, 
to drop to around 400 km s~^ again, at 5i?e- Thus, the 
rotation curve of this galaxy does not display the kind of 
"pseu do-Keplerian" fall-off of the that INapolitano et al] 
(|2007f ) identify with the "ordinary" class of ellipticals. 
Thus, our interpretation of NGC 4636 as a dark-matter 
rich galaxy is consistent with the latter class of galaxies. 
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